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The use of Lagrange multipliers in the context of quintessence/phantom scalar fields allows to con- 
strain the behavior of the scalar field, which provides a powerful tool, not only for the reconstruction 
of cosmological solutions but also for the study of some problems in cosmology and gravitational 
physics. In the present paper, we focus on the reconstruction of cosmological solutions capable of 
controlling the cosmological constant value by imposing a constraint on the scalar field, providing 
a relaxation mechanism of the value of the cosmological constant. The formalism is also extended 
to the study of phantom scalar fields with a future singularity and their conformal transformation 
to the Jordan frame, where a type of modified gravity, constrained by the Lagrange multiplier, is 
obtained. 

PACS numbers: 95.36.+X, 98.80.Cq 

I. INTRODUCTION 

In the last decade, a large number of Dark Energy candidates have been proposed in order to achieve explain 
the observational data, and also solve some old theoretical problems as the cosmological constant (cc) problem. 
However, the full understanding of the nature of Dark Energy is still very far of being complete. Firstly because 
of the impossibility to distinguish between the diverse theoretical models and the absence of new predictions, but 
also because of the complexity of some of these models, as modified gravity (for recent reviews on modified gravity 
and the reconstruction of cosmological solutions, see Nevertheless, it was proposed in a model which might 
provide the necessary tools for a better understanding of some of the existing problems in cosmological models and 
in gravitational physics. The model consists of a common quintessence/phantom scalar field, but as a novelty, the 
evolution of the scalar field has to satisfy a constraint introduced through a Lagrange multiplier, which basically 
behaves as an additional non-propagating scalar field, that retains a single dynamical degree of freedom. By a 
particular choice of the constraint, the differential equations for the two scalar fields, the quintessence/phantom 
field and the Lagrange multiplier, become first order equations, so that there are no propagating wave-like 
degrees of freedom. Such formalism can be easily extended to modified gravity instead of a scalar field, as it 
was proved in Ref. [3], which also provides interesting properties. This approach has shown that cosmological 
solutions which reproduce dark energy epoch can be easily reconstructed, but also it provides the necessary tools 
to solve the problems of a particular model, as violations at local scales or additional propagating modes of the 
graviton, which can be avoided by the constraint of the Lagrange multiplier (see 2]- 4]). Hence, the formalism can be 
used not only as a tool to reconstruct solutions but also to identify and isolate intrinsic problems of a particular model. 

The main purpose of this paper is the reconstruction of models capable to reproduce late-time acceleration, 
providing a natural explanation of why the value of the observed cosmological constant (cc) at the present time is so 
small compared with the vacuum energy density predicted by quantum field theories. By using the formalism of a 
scalar field constrained by a Lagrange multiplier, a relaxation mechanism of the value of the cc can be constructed, 
so that the cc problem can be dealt through the constraint on the scalar field, and a better understanding on its 
the evolution can be obtained (for other kind of relaxation mechanisms and their evolution, see ^] ) . Even, by such 
mechanism the inflationary epoch can be also reproduced, providing a unified picture of the entire Universe evolution. 
These results can be easily extended in terms purely of the metric tensor, as it was shown in Ref. Q, which may 
provide a more natural explanation, where the gravitational field is subject to some constraints. 

On the other hand, such model allows to study the behavior of some future singularities when a conformal 
transformation is applied, which relates the Einstein frame and the Jordan one. Specifically, we focus on the 
reconstruction of F{R) gravity from a phantom Universe described in the Einstein frame. Phantom cosmological 
models have been widely studied, as they carry very interesting physical information (see Refs. It is well 
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known that F{R) gravity is equivalent to a scalar-tensor theories class with a non-minimally coupling scalar field 
and null kinetic term, such that the action in the Einstein frame with a scalar field can be transformed into a f{R) 
action through a particular conformal transformation (see Refs. P-[loj). However, this is restricted to the cases of a 
quintessence field, free of future singularities (see Refs. [lH), while for the phantom case, the action becomes complex 
after being conformally transformed (see Ref. unless an additional component is included [l^. By using a 

Lagrange multiplier, a type of F{R) gravity can be obtained in the Jordan frame starting from a phantom solution 
in the Einstein frame. Hence, the conformal transformation of future singularities (for a classification of future 
singularities see Ref. p^). specifically the Big Rip singularity, is studied, which provides a better understanding of 
the future singularities when a particular conformal transformation is applied. 

The paper is organized as follows: in the next section, we introduce the quintessence/phantom models with a La- 
grange multiplier. Then, in section HI we reconstruct cosmological solutions with a relaxing mechanism for the cosmo- 
logical constant. In section IV, we reconstruct F{R) gravity in the Jordan frame starting from quintessence/phantom 
scalar fields in the Einstein frame. Finally, in the last section we analyze the transformation of a phantom field with 
a Big Rip singularity into /(i?) gravity in the Jordan frame. 



II. SCALAR-TENSOR THEORY WITH LAGRANGE MULTIPLIER 



In this section, we introduce the model where a minimally coupled scalar field is considered with the presence 
of a Lagrange multiplier that imposes a constraint on the scalar field. In the next sections, the reconstruction of 
cosmological solutions is done, where the Lagrange multiplier plays an important role. The general action considered 
here is given by, 

R Lu{(j)) 



S = / d^xV^ — - -^d.cj.d'^^ - V{(b) - \F{^, d.cj^d^^, ...) + Cm} ■ (1) 



2k 

Here = SttG and A is the Lagrange multiplier, while Cm refers to the matter lagrangian. The function 
F{(j),dfj,(j>d^(l), ...) is the constraint on the scalar field (j), while uj{(j>) and V{(j)) are the kinetic term and the scalar 
potential respectively. Note that the Lagrange multiplier behaves as a second scalar field, which does not propagate. 
By varying the action ^ with respect the metric g^^, the field equations are obtained. 



R^u 2 



where, 



tM = . (2) 

While the variation over A gives the constraint equation on the scalar field 0, 

^^(0, 9^09^0, ...) = . (3) 

Hence, solving the equation ([3]), the evolution of the scalar field (j) is fixed, which affects the Einstein field equations 
For some specific functions F, one can obtain interesting properties and cosmological solutions, where a relaxation 
mechanism of the vacuum energy density can be implemented, as shown below. In this paper, we are interested in 
studying spatially flat FLRW, whose metric is given by, 

3 

1=1 

Then, by introducing the metric ^ in the field equations ^ and specifying the function F(0, d^(t>d^(j), ...), the usual 
FLRW equations are obtained. Let us consider a simple example in order to show the effects of the constraint equation 



(5) 
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By solving this equation, the solution gives a constant scalar field (j>{t) = (f>o, the roots of the equation ([5]). Hence, 

the energy- momentum tensor gives an effective cosmological constant T^it'^^ = —5^1/^0, which can compensate for the 
large contribution of the vacuum energy density. Note that in general, equation ([5]) contains more than one positive 
root 00, so that each point provides an effective cosmological constant which may reproduce, or at least contribute, to 
the periods of accelerating expansion of the Universe history, unifying inflation and late-time acceleration under the 
same mechanism. However, here we are interested on those cases with a non-constant scalar field as more interesting 
results may be obtained, so from now on we consider the constraint equation. 



(6) 



As we are considering an homogeneous and isotropic metric such as FLRW, the scalar field is taken only as a function 
of time dependent </) = </'(t), so that the constraint ([6]) on the scalar field is. 



1 



= ±dt . 



2 vm^ 

By the metric (|4]), and using the condition ([7]), the FLRW equations are given by. 



(7) 
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(w(<;^) + 2A)C/(0) + V^(0)+^p, 



-2H ^ 



;(0)C/(0)-n'^)+^ 



(8) 



Here, we are considering matter species described by perfect fluids, where the subscript i refers to every kind of 
matter in the Universe with an equation of state given by pi — WiPi. Then, by using the above equations, the 
corresponding cosmological solution is obtained. The evolution of every matter specie can be easily calculated by 
solving the continuity equation. 



Pi 



3H{l + w,)p,=Q 



(9) 



where poi is an integration constant. Then, we can give a reconstruction method, similar to the case of free scalar 
fields (Ref. 0), in order to reconstruct the corresponding theory for some particular solutions H{t) and As the 
expression for (f) can be inverted, the time coordinate is written as a function of the scalar field t = t{(f>), and by using 
the constraint equation ([7]) and the FLRW equations ([H]), the theory that reproduce a specific solution is described 
by, 



C/(0) 



K 



1 r 



iH{t{4>)f + 2H{t{4>)) \ + ujU{(t>) + w.po^e 



^3(l+u),) /'<*' dt'H(t') 



(10) 



while the Lagrange multiplier yields, 



2ujU{(t){t)) + —H{t) + ^(1 + w.)po«e- 



'3(l+Wi) /* dt'H[t') 



(11) 



Note that uj{<j)) in the action ([T]) remains arbitrary, and we have assumed here a constant kinetic term a;((/)) = uj. 
Hence, we have implemented a method to reconstruct the corresponding action ([1]) for any given solution. Let us 
consider the example described by the solutions. 



c^[t) = , H{t) ^Ho + 



Hi 



(12) 



where Hq, Hi, and (po are constants. This solution describes a Universe that crosses through out two phases, the first 
one dominated by the second term in r.h.s. of the Hubble expression (fT2| . which produces a phase of decelerating 
expansion, identified with the radiation/matter dominated epochs. Then, a second phase dominated by the constant 
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-ffo, which describes a de Sitter solution at late times is reproduced. By using the expressions (|10l) . the corresponding 
potentials for the scalar field are obtained, 



Ui<t>) 



1 



2 ' 2 

While the Lagrange multiplier pT|) for this particular case is given by, 



-3Hiil+Wi) 



-3(l + i)i)-Ho ^ 



1 



Mi) =--2 



1 



i2 



(13) 



(14) 



Hence, we have reconstructed a theory, described by the action ([T]) and the potentials in ([T3|) . which reproduces the 
solution ((T^ . The behavior of the Lagrange multiplier is also fixed by (fHl) . Note that any kind of matter specie can 
be included in the equations (jl3p . We might include the vacuum energy density, whose EoS is given by the same one 
of a cosmological constant, i.e. Pvac = —Pvac, and its contribution to the dark energy density would be compensated 
by the presence of the scalar field. 

By this simple reconstruction method, we have shown that the constraint imposed on the evolution of the scalar field 
([7]) affects the expansion history and by a suitable choice of the potentials (ITOl) , the large expected value of the vacuum 
energy density can be relaxed by a natural mechanism. The next section is devoted to the analysis of some examples 
in more detail, where also the unification of inflation and dark energy is considered. 



III. RELAXING THE COSMOLOGICAL CONSTANT 



We are interested here to study the relaxation mechanism of the vacuum energy density in detail with explicit 
examples that are able to reproduce the dark energy epoch, and even the inflationary era. In the previous section, a 
general reconstruction method of cosmological solutions was presented, and some examples were studied. Nevertheless, 
here we are more interested in studying the evolution of the vacuum energy density, and the scalar held contribution 
to the dark energy density, and how the implementation of a relaxation mechanism can be realized by using the 
Lagrange multiplier. Let us consider the FLRW Universe described by the second equation in 

- ^(3i/' + 2H) = W,nP,n - Pvac + ^<P^ " V {(j)) , (15) 

here we have used the constraint equation Q, C/(0) = \(^. We have also included the energy density for matter p™, 
the vacuum energy density, p-uac — and we have assumed a constant kinetic term, given by a;(0) = 1. We define 
an effective cosmological constant in the equation (1151) that depends on time, and which can be written as. 

By the appropriate scalar potentials, we can implement a relaxation mechanism of the vacuum energy density. Then, 
the FLRW equation takes the form, 

- 31/2 _2H = n'^w.mPra - ^ef f , (17) 

Hence, the terms depending on (j) in (|16|) may act by relaxing the cosmological constant A, such that the effective 
value coincides with the one observed, Af,ff(to) ^ Hq, where the subscript o refers to the present epoch. In this sense, 
we will give some examples, starting with simple models where there is only a constant contribution to the effective 
value of the cosmological constant, to more complex models that are also capable of unifying the inflationary and 
dark energy eras. 



A. Example 1 



Let us consider a model where the potential U{(j3) is given by. 



[7(0)^^ + 1/(0), (18) 
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where is a positive constant. Then, the effective cosmological constant (fT6|) yields, 

Ae//=A-A^. (19) 

Hence, the value of the cosmological constant is compensated by the contribution of the scalar field. And in this case, 
the equation for the Hubble parameter (IT71) turns out to be the same as for the ACDM model, whose solution can be 
easily calculated. 

While the evolution of the scalar field can be calculated by introducing the potential ([T^ in the equation Q , and by 
specifying the potential V{(l)). If we consider for this case, V{(j)) = ^m^cj)^ with m being a constant, the solution for 
the scalar field is given by. 

On the other hand, the expression for the Lagrange multiplier can be easily calculated by introducing the above 
quantities in the first FLRW equation in ([8]), which yields A(i). Hence, this simple model shows a way to relax 
the value of the vacuum energy density by means of the scalar field, and the corresponding Lagrange multiplier. 
The evolution is the same as in the ACDM model, and the effective cosmological constant can be adjusted with the 
observational data, setting the value of A^. In the next examples, we study more complex models, where the effective 
cosmological constant can change with time, and the unification of the inflationary and dark energy epochs can be 
realized. 



B. Example 2 

We are now interested to consider a model where the effective cosmological constant (fT6|) is time-dependent. Let 
us consider now the potentials, 

f/(^) = 2F(0) = im2</)2 , (22) 

where m is a constant. It is straightforward to calculate the evolution of the scalar field by solving the constraint 
equation ([7]), which gives, 

^{t) = 0oe-^* . (23) 
Here = 8m^. Then, for this model, the effective cosmological constant ([TB)) is now given by, 

Ae// = A-A^e-"* , (24) 

where A^ = '^^'^ 0q, and we assume A^ < A. In this case, the effective cosmological constant evolves with time, 
whose behavior will depend on the value of a that determines the growth rate of the effective cc. In this model, the 
general solution of the equation (|17p is given by a much more complex function for the Hubble parameter than in 
the model analyzed above. However, for short time intervals, the solution can be well approximated by the ACDM 
solution, given in (pn|) . Note that the minimum value of the cc is given at t = 0, and the radiation/matter 
epochs can be reproduced. Then, the cc constant grows and the Universe enters in the dark energy epoch, as one 
expects. Note also that the value of a has to be very small in units of Gyrs~^ in order to match the observational 
value PK^fj ^ lO^*^^ GeV^ at the current epoch. 



C. Example 3 



As a last example, we want to consider a model capable to unify inflation and dark energy epoch under the same 
mechanism, assuming that the vacuum energy density dominates during a short time of period right after the origin 
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of the Universe providing a large value of the cosmological constant, then it decreases down a minimum limit until 
the present epoch, where it increases again. Hence, the U and V potentials considered here are, 

U{cb) = -, V{cf>) = 2.2(«02_^2) • (25) 

Here A^, 0O; a and to are constants. Then, by the equation ([7]), the evolution of the scalar field is given by, 

(j>{t) = 0oi ■ (26) 
While the effective cosmological constant is now described by, 

Ae/Z^A-A,!^. (27) 

This expression describes a cosmological constant that has a maximum value at t — Iq, where inflation occurs caused 
by the vacuum energy density, then the value of the effective cc decreases to a minimum, after which it starts to 
grow again. Hence, this model is able to reproduce the unification of the accelerating epochs of the Universe. It is 

straightforward to show that the minimum of the function (|27|) is given at tmin = to + y^tg + ^ , so that if we assume 

A^ ~ 2Aa (^to + \Jto ~^ J tti^ minimum value of the effective cc (p7)) can be fitted around zero at the present age 
of the Universe matching the small value given by the observational data. Following this example, the value of the 
effective cc should increase in the future, the time it takes to observe the possible effects depends entirely on a. Hence, 
an example that provides a relaxation mechanism for the value of the vacuum energy density at the present time but 
which also reproduces the early accelerated epoch, has been reconstructed. The results achieved here can be easily 
extended to modified gravities instead of the use of a scalar field, which can give a natural answer to the evolution of 
the Universe. 



IV. RECONSTRUCTING MODIFIED GRAVITY WITH LAGRANGE MULTIPLIERS 



We have shown that the use of Lagrange multipliers can provide a useful reconstruction technique, where interesting 
properties on the cosmological paradigm can be obtained. Particularly, the relaxation of the value of the cosmological 
constant can be easily achieved by using a scalar field constrained by a Lagrange multiplier. The same technique can 
be extended to modified gravity, as it was shown in Ref. 

In this section we are interested to transform the action ([T]) into the Jordan frame in order to obtain the corresponding 
F{R) gravity. It is well known that F{R) gravity is equivalent to a type of Brans-Dicke theory with a potential and a 
null kinetic term written in the so-called Jordan frame, such that F{R) gravity can be obtained from quintessence-like 
scalar fields by applying a particular conformal transformation on the Einstein frame (see Refs. [l0|-p^). Here we are 
interested to extend this conformal transformation to the case of a minimally coupled scalar field in the presence of 
a Lagrange multiplier. By applying a particular conformal transformation on the action ([1]), the corresponding F{R) 
action is obtained jTo| . 



r 

g^i, — il^gfj,i, where fl^ — exp ±k / < 



d(l)\ -uji4>) 



(28) 



Here the "tilde" on the metric refers to the Jordan frame. Note that the transformation ([28|) suppresses the kinetic 
term proportional to uj{(j)) in order to ensure that the corresponding Brans-Dicke-like action in the Jordan frame is 
equivalent to F{R) gravity (see Refs. @, @). However, the terms in front of the Lagrange multiplier remains also in 
the Jordan frame, where the action in vacuum is transformed as. 



Sj= I dx 



2k2 



R-V{^) - A 



1^2(0) 



(29) 



where V{(t>) = r2'*((^)U(0) and [/((/>) = il*((^)t7(0). By varying the action with respect to A, the constraint 
equation is obtained, similarly to the case of the above sections. Note that in general, the action is not equivalent 
to any F{R) gravity unless one imposes an additional condition on the conformal transformation (j28l) . By varying 
the action (l29t with respect to the scalar field, the equation for is obtained. 



^(0) 



dn{4)) R du(0) 



+ A 



(30) 
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Then, in order to obtain F{R) gravity in the Jordan frame, the last term in (l30l) have to become null, which can be 
achieved by imposing A = that eliminates the constraint on the scalar field and the action in the Einstein frame 
^ reduces to the usual action for a quintessence/phantom field, or we can impose that terms inside the brackets are 
null, which implies the condition. 

This condition basically restricts the form of the kinetic term uj{(j)), which remained arbitrary in the Einstein frame 
as showed in the previous sections. Here, it will be shown that F(R) gravity can be obtained from phantom universes 
in the Einstein frame by applying the conformal transformation obtained, contrary to the usual phantom case (see 
p^). Hence, by the condition (l3Tt . the equation (|30)) yields. 



^(^)^4-^-0 ^ c^^m (32) 



By solving this equation, we can obtain the corresponding action for F{R), 



Sj — I dx^ \/ —I 

where. 



F{R) + A {^^^^^^R^i^R + V{R) 



(33) 



F{R) = ^^R - , and W[R) ^ n^m) ' . (34) 

The action (1551) corresponds to a type of F{R) gravity introduced first in ^3], where the scalar Ricci 7?, and for 
instance the metric tensor, is constrained by the Lagrange multiplier, similar to the case of scalar fields studied above. 
It is straightforward to show that assuming a flat FLRW Universe in the Jordan frame, a relaxation mechanism of 
the cc can be reconstructed starting from the action ([M)) . but in this case, the relaxation is realized in terms of the 
metric instead of scalar fields, what can provide a more natural explanation to the observed value of the cosmological 
constant. However, here we are interested to show the mathematical relation between the action defined in the Einstein 
frame in terms of a scalar field ([T]), and the F{R) action obtained in the Jordan frame (fSS)) . and specially the relation 
between the corresponding cosmological solutions. Note that we are referring here to the mathematical equivalence 
between both frames which may help to understand their physical relation, however, the physical equivalence of both 
frames is a wide open discussion beyond the aim of this paper (for a review see Ref. [l5|). It is well known that 
phantom scalar fields defined in the Einstein frame has not a F{R) correspondence in the Jordan frame (see Ref. [T^). 
unless additional components are assumed (see Ref. ^T3\). However, in the case of the action ([T]), the conformally 
transform to F{R) gravity is not restricted only to the non-phantom case. Let us first write the relation between the 
cosmological solutions for each frame. By assuming a flat FLRW metric in each frame, the metric components and 
coordinates in both frames are related by means of the conformal transformation ()28|) , 



a(t) ~ dt , 

Let us consider now an explicit example, where the corresponding F{R) action is obtained as well as the cosmological 
solutions in the Jordan frame starting from the Hubble parameter in the Einstein frame, 

H = H{t) , 0(0 = cj^ot . (36) 

,2 

By the previous section the solution for the scalar field ([55)1 corresponds to a constant potential U {(p) — ^ , according 
to the constraint equation ([7]). We can now rewrite the equation pip in terms of the variables defined in the Einstein 
frame according to the transformations (1351) . 



nW4>-^'P^ + 3nH<p + AU{c^)^ + n^^^o . (37) 

dip dip dip 

Here the dots on the scalar field refers to derivatives with respect the time t defined in the Einstein frame. By 
assuming the solution (1551) . the equation (1571) is considerably simplified, 

</)o(2-0o)5r + 30oiff^(</') = . (38) 
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By conformal transformation defined in ([28]) . the equation (|38|) can be solved for the kinetic term Ld{(j>), which gives 

27 

8k^(2 - 0o) 

Hence, by introducing this function in the conformal transformation (j28p . we are able to obtain the corresponding 
F{R) gravity in the Jordan frame Note also that for this explicit example, uj{(j)) > independently of the 

cosmological solution H{t) in the Einstein frame, such that a phantom Universe can be conformally transformed in 
terms of a F(R) action in the Jordan frame, which supposes a difference with respect to the usual phantom scalar 
fields |12|]. 



V. PHANTOM FIELDS AND FUTURE SINGULARITY IN THE CONFORMAL FRAME 

Let us now analyze the transformation of an explicit Hubble parameter that contains a Big Rip singularity in the 
Einstein frame. 



H{t) 



t.-t 



a(t) 



ao 



t, - t 



(40) 



where oq is an integration constant and is the Rip time. The solution ()40|) describes a super accelerating Universe 
in the Einstein frame that ends in the so-called Big Rip singularity at t — ts- We are now interested to see how the 
singularity is transformed in the Jordan frame ((33)) . By introducing the expressions ([36]) and (|40)) in ([39)) . the kinetic 
term takes the form. 



27 



8K2[(2-0o)(t.-0/</'o)]2' 
While the scalar potential V{(j)) for the solution (HUj) can be obtained by using the expression pi7)) . 

Vo 



(41) 



(42) 



where Vb = ^^'^x^^l^y'^^^ ■ The conformal transformation ([25)1 is given by. 



f]2(0) 



Hence, the cosmological solution fj40|) is transformed in the Jordan frame through ([35]l. which yields 



a dt t 



(43) 



(44) 



where gq — ^ ■ This solution describes an expanding or contracting Universe depending on 

the sign in (j4T|) and on the value of the constant (pQ. Nevertheless, independently of the type of the cosmological 
evolution, the solution pi)) contains a singularity at t = that can be identified as a Big Bang (Crunch). Therefore, the 
singularity in the Einstein frame (the Big Rip singularity) is transformed in the Jordan frame as an initial singularity, 
which means that the conformal transformation (|28p does not remove the singular point but displaces it. 
In order to reconstruct the F{R) action in the Jordan frame, one has to solve the equation (|32p. For the example 
studied here, the solution is given by. 



-1+ 



4>{R) = (f>0 ts - (ai?) ±3*0+2(2^*0) 



And the F{R) action in the Jordan frame (155)) yields. 



where a = ± 



3</>o 



(45) 



F{R) 



2(.»0-2) 



(aR) ±3*0+4-2*0 

2^ 



Voa 



R , W{R) = !3{aR) ±3*0+4-2*0 



(46) 
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where P = ig^4ya(-|-0p^^4^2<^o)^^('±0o+2-0o)^ ' Hence, we have reconstructed the corresponding F{R) gravity in the 
Jordan frame, whose cosmological sohition is given in (|44p . Therefore, an action with a phantom scalar field constrained 
by a Lagrange multiplier in the Einstein frame ([B]) can be conformally transformed to a F{R) action in the Jordan 
frame. It is important to point out that the singularity remains in the transformed frame, where appears as an initial 
singularity. 

VI. DISCUSSIONS 

In the present paper, we have studied the use of a Lagrange multiplier to reconstruct cosmological solutions, and 
to study some problems related to gravitational physics. We have shown that the usual quintessence/phantom scalar 
fields, restricted by the presence of a Lagrange multiplier, can reproduce any kind of cosmological solution, while the 
evolution of the scalar field can be maintained under control, and useful information on the cosmological evolution 
can be analyzed. In this context, we have proposed a relaxation mechanism of the cosmological constant (cc), which 
provide a way to study the evolution of the cc value. Hence, a scalar field with no additional propagating wave modes, 
and whose equations are simpler than the usual case can parametrize and compensate the value of the vacuum energy 
density in a natural way, providing a method for a better understanding of the cc problem, whose cosmological 
evolution can be easily studied. In this sense, we have provided several examples in order to show the evolution of the 
"dark energy" density, and how the value of the vacuum energy density is reduced. The next step will be a deeper 
study of this kind of models, and its constraints with the observational data, but this is beyond of the purpose of this 
paper. In addition, we have shown that solutions capable to unify inflation and late-time acceleration under the same 
mechanism can be reconstructed. Naturally this approach is extended to F{R) gravity, where the restriction affects 
to the metric tensor. Such constraint is expected to become important at large scales, while at small scales the usual 
equations of General Relativity should be recovered. As in the case of a scalar field, the Lagrange multiplier is a pow- 
erful tool to control the problems of a given model, and it gives important information on the evolution of the solutions. 

On the other hand, the formalism has been also used to study the behavior of future singularities, specifically the 
so-called Big Rip singularity, when a conformal transformation that led to a type of F(R) gravity in the Jordan 
frame is applied. It is well known that the reconstruction of F{R) gravity is technically not possible in the usual case 
of phantom scalar fields. However, in this case, we can have a positive kinetic term of the scalar field defined in the 
Einstein frame due to the constraint equation provided by the Lagrange multiplier, which reproduce the phantom 
evolution in the original frame. Hence, this allows to study the behavior of future singularities when this kind of 
conformal transformation is applied. We have shown that the singularity in the Einstein frame remains also in the 
Jordan one, but with a different behavior. By allowing the time coordinate scale with the conformal transformation, 
we have shown that the Big Rip singularity is transformed into a Big Bang (Crunch) singularity. This may help for 
a better understanding of the physical relation between both frames, however such discussion is beyond the purpose 
of this paper (see Ref. [lB|)- Other aspects about the relation between both conformal frames and its (un)physical 
equivalence are recently studied in Ref. (l6j . 

Hence, the use of Lagrange multipliers provides a scalar field model with simpler equations and without additional 
wave modes, that can be used to study the evolution of the cc value and provide a way to analyze its evolution, 
and reconstruct the entire evolution of the Universe in a natural way. In addition, we have shown that F(R) gravity 
can be restored from a model containing a Big Rip singularity, which gives a better understanding on how future 
singularities behave. 
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